Giant vortices in small mesoscopic disks: an approximate description 
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We present an approximate description of the giant vortex state in a thin mesoscopic super- 
conducting disk within the phenomenological Ginzburg-Landau approach. Analytical asymptotic 
expressions for the energies of the states with fixed vorticity are obtained when a small magnetic 
flux is accumulated in the disk. The spectrum of the lowest Landau levels of such a disk is also 
discussed. 

PACS: 74.20 De, 74.60 Ec 



I. INTRODUCTION 



Progress in microfabrication technology has made mesoscopic superconductors (mesoscopic disks, in particular) 
a very popular subject of study (see, for example, Refs. Mesoscopic samples have sizes comparable to the 

coherence length ^ and the magnetic field penetration length A. The behaviour of such structures in an external 
magnetic field H is strongly influenced by the sample shape |^ and may lead to various superconducting states [^,|| . 

The vortex structure in mesoscopic disks with finite thickness was studied early using a numerical solution of 
the system of two coupled non-linear Ginzburg-Landau (GL) equations p. For thin disks it is possible to simplify 
the solution of the problem by averaging the order parameter over the disk thickness and by representing the 
superconducting order parameter as a superposition of eigenfunctions of the linearized GL equation 1^,^,0. But, even 
then the problem has still to be solved numerically because of the presence of confluent hypergeometric functions. 

It was shown previously that in some limiting cases the description of the vortex structure can be simplified. 
Buzdin and Brisson described certain vortex configurations in a small superconducting disk within the London 
approximation. Akkermans et al. calculated the vortex structure of mesoscopic disks at the vicinity of the dual point 

= 1/2 and in the London limit [|o|. Recently, the structure of a giant vortex in an infinite plane was analyzed 
analytically within the GL theory for arbitrary values of k in the limits of small and large values of the vorticity Jni | . 

In the present paper we obtain analytical results for the free energy of small disks within the approach of Refs. [^^|, 
using the magnetic flux trapped by the disk as an expansion parameter. 



II. GENERAL THEORETICAL FORMALISM 



We consider a mesoscopic superconducting disk with radius R and thickness d << X,^ magnetized by the external 
magnetic field H = {0,0, H) which is uniform and directed normal to the disk plane. The theoretical model was 
already described in detail in Refs. and therefore we scetch only those steps which are necessary for our analytical 
approach. For a thin disk, to a first approximation, the magnetic field is uniform inside the disk and equal to the 
external one. As a result of this approximation, the distribution of the superconducting order parameter in the disk 
plane V'(p) is described by the first GL equation 

(-ZV2I3-1)'V^ = V^(1-|^|') , (1) 

with A{p) ~ (0, Hp/2, 0) and with the boundary condition at the sample surface 

(^iV2D -A)^ =0. (2) 

V / p=R 

The index 2D refers to the two-dimensional operator. Due to the circular symmetry of the sample we use cylindrical 
coordinates: p = {p,0) {p is the radial distance from the disk center, 9 is the azimuthal angle). All distances are 
measured in units of the coherence length ^, the magnetic field in Hc2 = 2^/^kHc, where He is the thermodynamical 
critical field. 

In the giant vortex state the order parameter can be expressed as 



^{p,e)=[-Kj-^] fL{p)e^p{iLe), (3) 
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are the eigenfunctions of the Uncarizcd Eq. dll) and 
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determines the spectrum of the lowest Landau levels, 
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is the Kummer function p^ . The value of v is determined by a non- linear equation, which results from the boundary 
condition (||) 



[L - $) M {-V, L + 1, $) - {-V + 1, L + 2, $) = 0, 

_L + 1 



(8) 



where $ = HR^ /2 (measured in units of $0 = TiTic/e) is the magnetic flux through the disk in the absence of any flux 
expulsion. The free energy, measured in = H'^V/8tt units, is 



F 



(9) 



III. APPROXIMATE RESULTS 

First, we have to solve the non- linear Eq. |8| In the case of small (f> we can restrict the series (^ for the Kummer 
hmction by the first four terms. The solution of Eq. ^ can be written as an infinite series of which the first four terms 
are 



{L + 1){L-<P) 
2$ 



1 



L + $ 



L (L + $) {L + $) ($2 _ L2 + 2L3 



2{L + 2) 2(L + 2f(L + 3) 4 (L + 2)' (L + 3) (L + 4) 
In the limit of large <I> >> 1 we use the asymptotic expression for the Kummer function 

M (a, c, z) = exp(z) [l + O (izp^)" . 

Substituting it in Eq. || we obtain the equation (L + ^)/T{—j^) = 0, which leads to the solution 

t'inf = 0. 



(10) 



(11) 



(12) 
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A. Lowest Landau levels 



Substituting v from Eq. (10) into Eq. (0) we obtain the eigenvalues A in the <i> << 1 limit: 
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(13) 
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Because series (|^ is asymptotic, each new term gives a contribution to all the coefficients of the <i>-series. 
In the <i> >> 1 limit we obtain for A the asymptotic expression 

Ainf = -1 + ;^, 



(14) 



which does not depend on the vorticity L. In Fig. |l] the A($) dependences are shown for a disk with radius R = 3^. 
Only the curves with L < 6 are shown. The solid curves represent the dependences calculated from the numerical 
solution of expressions (^) and (^) and the dashed curves are our new approximate results as given by Eq. (|l^). Notice 
that for small values of the magnetic flux the approximate expression (13) describes quite well the lowest Landau 
level spectrum. 




FIG. 1. The magnetic flux dependence of the lowest eigenvalues of the linearized GL equation (|l|) for different vorticities L. 
The "exact" numerical results are shown by the sohd curves. The approximate results are shown by the dashed curves ($ << 1 
limit) and by the dash-dotted curve ($ >> 1 limit). 



B. Free energy of giant vortices in the small $ limit 

The eigenvalues A, Eq. (^), determine the minimal free energy F, Eq. (|^), of the giant vortices. We expand Eq. (|^) 
in a series with respect to $. We found that for small disks with radius R — {2 ~ 3)^ it is sufficient to keep in the 
series (M) only the flrst four terms. Expanding the free energy up to order <I>^ we obtain the approximate expressions: 
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$2 / 557 1 1 \ 4 

^ 7.87 13.34 / 25.66 51.38\ ^ / „ 139.08 263.87\ 

96.56 - + ™^ - ^501.67 - ^ + '-^] 



R'^ J \ ' R"^ R^ 

The free energy for the next L states has the same form as -Fl=i but with different numeric coefficients and, therefore, 
will not be give here. 




FIG. 2. The free energy of the giant vortex states for a thin disk with radius 7?/^ — 2.0. The solid curves are the numerical 
solution and the dashed curves are from our analytical approximation. The dash-dotted curves are from our analytical expansion 
up to the $^ term. The dotted curves are from the analytical expression of Ref. . 

In Fig. H the free energy of the disk with i? = 2^ is shown. Notice that limiting oneselves to terms in F up to 
leads to a poor approximation, while including terms up to results in an excellent agreement with the numerical 
results. The dotted curves in Fig. ^ present the free energy calculated with the expression 



proposed in Ref. ||13[ (and rewritten in the notation of the present paper). Notice that expression ( |17| ) describes less 
accurately our numerical result (solid curves). It agrees for L — with our $2 approximation, but for L = 1 there is 
a substantial deviation with our approximate result. This disagreement is even larger for L = 2 and is therefore not 
given. 



IV. CONCLUSION 



We calculated approximate expressions for the lowest Landau levels and the free energy of the giant vortex states 
in thin mesoscopic disk with small radius. We found that these approximate analytical expressions agree very well 
with the more involved numerical calculation (EUtI and are a substantial improvement to those found in Ref. p3||. 
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